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Abstract 

We study the diffractive photo- and leptoproductions of J /ip and tp' on the 
proton, and examine the validity of the small-size color dipole approxima- 
tion to the production of radially excited heavy quarkonium. The dipole 
model predicts a small ratio of tp' to J/ip photoproduction cross sections, 
which does not agree with experimental data. We show that this discrepancy 
originates from a large transverse size of ip' which makes the convergence of 
the transverse size expansion questionable, and the calculation without the 
dipole approximation turns out to be consistent with the data. Productions 
of T'(2S, 35) are also discussed, and the dipole approximation is found to be 

reasonable for the T-family. 
13.85.Ni,14.40.Gx,12.38.Bx,12.39.Jh 

Typeset using REVT£X 



-mail : ksuzuki@nt . phys . s . u- tokyo .ac.jp 
-mail: arata@nt. phys. s.u-tokyo. ac.jp 
-mail: itakura@rcnp.osaka-u.ac.jp 
-mail : alam@nt . phys . s . u-tokyo .ac.jp 
e-mail: hatsuda@phys. s.u-tokyo. ac.jp 



1 



I. INTRODUCTION 



Recently diffractive phenomena in QCD have attracted considerable theoretical and ex- 
perimental interests. In particular, diffractive photo- and leptoproductions of vector mesons 
off the proton, 7M + p — > V(p°,uj,(j), J/if>, . . .) + p, provide crucial constraints on proper- 
ties of the QCD Pomeron and the vector meson wave functions [ffl-Rj. In this report, we 
focus on the heavy quarkonium production, where the heavy quark mass plays a role of the 
hard scale. Due to the Regge kinematics, the vector meson production can be understood 
as the scattering of the color dipole with the target nucleon. In this case we expect the 
amplitude to be dominated by the four diagrams illustrated in Fig.l. The upper part of 
each diagram, where the heavy quarks interact with two-gluons, is assumed to be calculated 
within perturbative QCD even for the photoproduction. On the other hand, the lower blob 
represents both non-perturbative structure (impact factor) of the proton and complicated 
gluon dynamics at small-x which is identified with the QCD Pomeron. 

The color dipole approach is adopted to evaluate the J /if) production in refs. [|2|||f|, and 
gives reasonable results in agreement with the data. However, as pointed out by Hoyer and 
Peigne very recently, the ratio of if)' to J/ if) photoproduction cross sections calculated in 
the dipole approximation is much smaller than the measured values [RJ at HERA, when one 
uses typical non-relativistic quark models to obtain the charmonium wave functions. The 
formula used in the dipole approach |2|||[| is obtained by assuming that both transverse 
momentum of exchanged gluons and the transverse size of the cc pair are small. However, 
the non-relativistic quark model indicates that the typical size of if)' ~ 0.8fm is twice as 
large as the radius of J/if>. Hence, the validity of the small-size dipole approximation is 
questionable for ip'(2S) and other mesons with the large radii. We may expect considerable 
contributions from the scattering of the color dipole with a large transverse fluctuation, if the 
7 and cc wave functions allow such a large separation between the quarks. Above argument 
motivates us to calculate the diffractive J /if) and if)' productions without resorting to the 
dipole approximation. 
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II. CALCULATION OF DIFFRACTIVE HEAVY QUARKONIUM PRODUCTION 



Let us briefly introduce the framework to calculate the diffractive heavy quarkonium 
production, which was developed in ref. 0. We consider the forward heavy vector meson 
production 7 < ■*' ) (q , ) + p(p) — > V(q) + p(p) shown in Fig.l with the kinematical condition 
W 2 = (p + q) 2 ^> Q 2 {= —q 2 )- In this small- x kinematics, the lifetime of a color-singlet 
quark- ant iquark fluctuation is large compared to the typical interaction time scale [il . 
Then the production amplitude can be expressed as a convolution of the gg-proton cross 
section with the wave functions of initial (7) and final (V) states. In terms of light-cone 
wave functions of the photon and vector meson ipv, the production amplitude at t = 
is proportional to 

A(z) = Ana s J d 2 b J dl 2 ^^- ^(Q 2 ; z, b)* [1 - J (l t b)} VvO, 6) , (1) 

where z and b are longitudinal momentum fraction of a quark and transverse distance 
between the cc pair, respectively, x is the longitudinal momentum fraction of the gluon 
x pa (Q 2 + 4m 2 ) /W 2 with m being the heavy quark mass, and U the gluon's transverse 
momentum. f(x,lt) denotes the unintegrated gluon number density of the proton, which 
contains a complicated non-perturbative gluon ladder. a s and Jo(bl t ) represent the QCD 
coupling constant and the ordinary Bessel function of the first kind. After multiplying suit- 
able functions of z and kinematical factors, and carrying out ^-integration, we obtain both 
longitudinal and transverse differential cross sections for the heavy quarkonium production 



' da^ 
~dt 



oc 

t=0 



M 2 



(2) 



where My is the mass of the heavy quarkonium. 

In eq. ([D the photon wave function %p 1 with 4-momentum square Q 2 is calculated pertur- 
batively, while ipv is a non-perturbative object. We shall use the non-relativistic potential 
model which successfully reproduces the properties of the heavy quarkonium to evaluate ipy 
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In the following calculations, we use two types of the potential; the QCD motivated poten- 
tial which essentially consists of linear and Coulomb potentials ||12|| , and a phenomenological 
logarithmic potential |13| . Our procedure to evaluate the light-cone wave functions is in line 



with those of refs. |||9|]. Solving the Schrodinger equation with a potential V(r), we obtain 
non-relativistic wave functions ipy{r), and also subsequently those in the momentum space 
ipv{k) through the Fourier transformation. We make a simple kinematical identification 
between light-cone variables (z, k t ) and k = (k z , k t ) as 



1+ h * 1 , (3) 

k 2 + m 2 



1 

z — - 
2 

where ±k t are the transverse momenta of the quarks. This relation gives the light-cone wave 
function; 



/ k 2 + m 2 \ 1/4 
ijj v (z,k t ) = [—^- — -ipvWkl 



k 2 t + (2z - l) 2 m 2 



(4) 



A[z(l-z)} 3 ) rv y 1 \ 4z(l-z) 

Finally, performing the two-dimensional Fourier transformation in the kt space, we obtain 
the light-cone wave function ijjy(z,b). 

If we take the infinitely heavy quark mass limit, where the Fermi motion of heavy quarks 
in the quarkonium can be entirely neglected, i.e. ipy — ${ z ~ 1/2)5^(6), the production 
cross section is simply related to the leptonic decay width of the heavy quarkonium [|TJ. 
However, for the charmonium, there exist substantial Fermi motion corrections at small Q 2 



Assuming that (bl t ) <^ 1 and keeping terms up to second order of (bl t ) PUd; one can 
rewrite eq. (0) as 



A(z) = ot s xG(x, Q 



eff) 



7T / d 2 b^(Q 2 ;z,byb 2 ^ v (z,b) 



(5) 



where Qlff = z(l — z)Q 2 + m 2 . In eq. (|) the unintegrated gluon density has been identified 
with the gluon distribution function measured in DIS via 

, dxG(x, I 2 ) 



f(x,l t ) = l 



1 dl 2 



The approximation (|5|) is often called the color dipole approximation. It gives a simple 
physical interpretation of this process, in which the scattering amplitude of the qq dipole 
with the target has a geometrical expression, a q g oc b 2 . In addition, the factorized formula 
eq. (|5|) makes it possible to constrain the gluon distribution function of the proton from 
the diffractive vector meson production. Typical results for J/if> production cross section 
obtained by using eq. (f|) are available in refs. [|3lJl4|, and in agreement with the data. 
However, if we focus on if}', which is a radially excited 25"-state of the charmonium, the 
calculated ratio of if>'(2S) to J/ip(lS) cross sections at the photoproduction point Q 2 = is 
much smaller than the experimental data as pointed out in ||. The ratio calculated using 
eq. (^|) is shown by the dashed curve in Fig. 2. Although the theoretical prediction agrees 
with the data at larger Q 2 , it has a factor 3 discrepancy at Q 2 = 0. We have tried to 
evaluate it with several potential models, but none of them can describe the data as far as 



the parameters of the potential are fixed to reproduce the charmonium properties [14 



Let us now examine carefully the dipole approximation which leads us to eq. ([5]). For the 
leptoproduction with large Q 2 eq. ([5]) works well, since the photon wave function projects 
out the cc fluctuation with the small transverse size. This approximation is also reasonable 
for the photoproduction of the IS" heavy quarkonium due to the small transverse size of the 
J /if}. However, for the photoproduction of the radially excited vector mesons, it is not clear 
whether the small-(6/ f ) expansion of eq. (||) is valid. Since the size of ip'(2S) turns out to be 
about 0.8fm, twice as large as that of J/ip(lS), in the typical non-relativistic quark model, 
the color dipole expansion may not work in this case. Namely, there may exist non-negligible 
contributions to the integral eq. (|l|) ^from an overlap of the large-size cc configuration and 
the if)' wave function. 

In this work we explicitly evaluate eq. ([[]) instead of using the simplified formula (|5|). To 
do this, we need the explicit form of the unintegrated gluon density f(x,l t ), which is not 



precisely known. Here, we follow the prescription of refs. [jnj,f|] to carry out the /(-integration 
in eq. ([!]) by relating f(x, l t ) to the derivative of the gluon distribution; 
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JO 



k J o h J Ql k 

Q'tff 2 dxG{x,lt 



xG(x, Ql)A(l t = 0) + / J" dl f X ^ lt) A(l t ) (6) 



where A{l t ) is a given function of l t in eq. (]l|), and Ql plays a role of the infrared separation 
scale. 

III. RESULTS AND DISCUSSIONS 

Using the CTEQ4LQ |16[ parametrization for the gluon distribution and taking Qq = 
0.8GeV 2 , we obtain the result shown in Fig.2. Calculation without (with) the dipole approx- 
imation is depicted by the solid (dashed) curve. For large Q 2 the difference between them 
becomes small as expected. However, around Q 2 = 0, the calculation of eq. ([!]) with eq.(^) 
clearly larger than the result obtained by eq. (^) by factor about 2 |T7[ , and shows a reason- 



able agreement with the experimental data ||10| . We demonstrate the major improvement 
of the ratio is caused by the contribution from the large-size qq fluctuation. There may be 
several other mechanisms to enhance the ip' to J/ip ratio further, which will be discussed 



elsewhere 14 



If we vary the value of Ql from 0.4GeV 2 to lGeV 2 , resulting ratio changes within a few %. 
Our results are also found to be insensitive to the choice of the potential model. We show in 
Fig. 3 results with GRV98LO [|18j parametrization for the gluon distribution function (thin- 
solid curve) in comparison with those with the CTEQ4LQ distribution (thick-solid curve). 
Calculated ratio is slightly modified due to the difference of the gluon distributions between 
CTEQ4 and GRV98 at low Q 2 . 

For completeness, sensitivity of the results to Q 2 // i s demonstrated in Fig. 3, because 
several choices of Qlff are adopted in the literature fHHH. In Fig. 3 results with Q 2 ff — > oo 
as an extreme case are shown by the thick-dashed (CTEQ) and thin-dashed (GRV) curves, 
while those with Qlff = Q 2 /4: + m 2 by the solid curves. Although increasing the upper limit 
of the l t - integration improves the result, significant enhancement at Q 2 = comes from 



6 



higher order terms of b 2 in eq. (|l|). 

We also consider the production of T(1S, 2S, 3S')-states. One can naively expect that 
the transverse separation b is small enough to ensure the validity of the dipole expansion 
for the bottom quark case. Results for T'(2S)/T(1S) and T'(3S)/T(1S) are shown in 
Fig. 4. They are consistent with those obtained in ref. |T9| by the dipole formula. There are 



only slight differences between the calculations without (solid) and with (dashed) the dipole 
approximation even for the T'(3S) case, whose radius is about 0.6fm. 

In conclusion, we have studied the diffractive photo- and leptoproduction of radially 
excited heavy quarkonia, and discussed the validity of the small-size color dipole approxi- 
mation, which gives the smaller value of the ip' to J / ip ratio compared with the experimental 
data. It is pointed out that the contribution from the higher order terms of b 2 in the qq- 
nucleon scattering amplitude, which is inherent in the original formula eq. ([I]), is important 
for the radially excited ip' case. We have demonstrated that the calculation of eq. (|l|) without 
the dipole approximation provides a much larger value for the elastic photoproduction ratio 
da(ip'(2S))/da(J/i(j(lS)) than the dipole result, and is consistent with the HERA data. As 
claimed in ref. ||, relativistic treatment of the charmonium system may be of some im- 
portance to construct the ip' wave function. However, what we have shown here is that, 
even within our non-relativistic framework where there are some ambiguities of choosing the 
kinematics e.g. in eq. @, substantial improvement for the ip' production can be achieved 
by evaluating eq. (JJ) instead of eq. Our results suggest that a careful treatment is nec- 
essary to study the radially excited light vector meson leptoproduction, e.g. p'(2S), <fi', ■ ■ ■ at 
moderate Q 2 . We have also discussed the T'(2S, 3S) production processes, and found that 
eq. (HD is a good approximation for the bb system because of its heavier quark mass. 
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FIG. 2. Ratio of ip' to J/ip production cross sections at W = 90GeV. Results obtained by the 
calculations with Qq = 0.8GeV 2 and the color dipole approximation are depicted by the solid and 
dashed curves respectively. 
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FIG. 3. G(x,Q 2 ) and Qgjj-dependence of the results with Qq = 0.8GeV 2 . See text for details. 



0.8 



0.6 



0.4 



0.2 



0.0 







2S full 
3S full 
2S Dipole 
3S Dipole 



10 



20 



30 



Q 2 [GeV 2 ] 



40 



50 



FIG. 4. Ratios of T'(2S, 3S) to T cross sections. Thick solid and dashed curves denote 
da(T'(2S))/da(T(lS)) obtained by the full calculation and dipole approximation respectively. 
Thin solid and dashed curves indicate the corresponding results for T'(35)/T(15). 
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